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EARTHQUAKE STRESSES SHEAR BUILDINGS 


The elastic stresses multistory building due earthquake displace- 
ments its foundation are determined means the normal modes 
vibration the building and numerical integration. The building assumed 
move due the deflection its columns, and the mass the structure 
concentrated the floors, which are infinitely rigid. The rocking the build- 
ing elastic soil and the influence internal damping are taken into ac- 
count. The stresses the multistory building are compared with those 
obtained assuming the building act cantilever shear beam. 
numerical example illustrates the application the method and permits 
evaluation its accuracy and efficiency comparison with rigorous solution 
the same problem. 


INTRODUCTION 


The dynamic behavior elastic framed building acted upon impressed 
displacements (earthquakes) time-varying forces (blasts) great practical 
importance. Unfortunately, the mathematical difficulties involved the 
solution dynamic problems have prevented rational design such struc- 
tures; hence, any step the simplification their solution 

The dynamical analysis framed structures considerably simplified 
whenever the structure has height-to-width ratio and whenever its 
columns are relatively flexible, which case the structure will called 
“shear shear building possible neglect the influence 
direct stresses the columns and assume that: (1) The motion the build- 
ing due the bending deflection the columns only; and (2) the mass 
the building lumped the floor levels. Within the limits elastic behavior, 

comments are invited for publication; the last discussion should submitted 
September 1953. 

Associate Prof. Civ. Eng., Columbia Univ., New York, Y., and Cons. Engr., Ammann 


Whitney, New York, 
Bracing,” Spurr, McGraw-Hill Book Co., Inc., New York, Y., 1930. 
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the building acts elastic system with many degrees freedom there 
are floors, and the determination the stresses reduced the integration 
set simultaneous ordinary differential equations. 

what follows, method for the solution this problem presented for 
the case earthquake stresses. The same procedures could applied, 
practically without change, the blast problem. 

Notation.—Letter symbols this paper are defined where they first appear, 
the text diagram. the symbol used subscript, re- 
fers the main symbol its appropriate floor; and the symbol used sub- 
script, denotes one the vibrational modes the structure. Single 
double dots over any symbol, denote differentiation with respect time. 


2 
= 
a 
© 


Consider slice framed building, one unit wide, with floors. Fig. 
let: 
the total horizontal shear force the ith floor; 
the story height between the 1)th and the ith floors; 

the spring constant under the ith floor—that is, the shear the ith 
floor necessary displace that floor horizontally one unit 
relative the 1)th floor; 

the modal displacement the ith floor; 
the absolute horizontal displacement the ith floor time and 
the impressed horizontal motion the foundation the building, 


considered identical with the earthquake displacement compon- 
ent the plane the frame. 


-- 
Fre. 1 
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The equation motion the ith floor obtained from Newton’s second 
law ma) the form: 


which the dots connote differentiation with respect time. 

The free vibrations the structure are characterized the fact that each 
floor moves simple harmonic motion. The structure can vibrate freely 
different modes defined the equation 


which X(i) the displacement the ith floor the mode frequency 
Substituting Eqs. Eqs. the ordinates are found satisfy the system 
homogeneous linear equations— 


which 


Eqs. have roots different from zero only their determinant vanishes, 
thus: 


quencies the modes. Once the values have been found, Eqs. deter- 


mine the the jth mode within multiplying constant 
This constant may chosen make 


N N 
i=] 


im] 


which case the modes 


Methods Th. von and Biot, McGraw-Hill Book 
Co., New York, Y., 1940, Chap. 


and 
— 
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—are said normalized. can that the modes (or are 
orthogonal; that is, that 


The displacements the structure due any impressed displacement, 
applied force, can expressed? series products the modes and 
unknown functions called the “generalized co-ordinates” the struc- 
ture. the study earthquake displacements convenient expand, 
terms generalized co-ordinates, the displacements the ith floor relative 
the foundation—that is, set— 


j= 
substitution Eqs. the expressions for the kinetic energy: 


the structure, and Eqs. and may seen that: 


N N 
= 


N 
which the total mass the building and 


N 


The equations motion the structure terms its generalized co- 
ordinates are Lagrange’s 


where the “generalized given, terms the force applied the 
ith floor, by: 


and 
and 
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Eq. could used determine the force exerted the ground the 
foundation, but will not considered this paper. the present case the 


only external force applied the foundation and but not 


function the and hence all the are zero. means Eqs. and 10, 
the Lagrangian equations become 


The acceleration graphs supplied strong-motion records‘ are jagged 
that often inconvenient integrate Eqs. 14. Then becomes practical 
introduce new generalized co-ordinates defined the equation: 


terms which the equations motion become 


and contain the displacement zo, forcing function, rather than the ac- 
celeration Eqs. and demonstrate the advantage gained the use 
generalized co-ordinates and normal modes. The burdensome integration 
the original system simultaneous equations (Eqs. has been reduced 
the much simpler task integrating the separate Eqs. 
over, will often found sufficient from engineering point view take 
into account only the first few (two three) modes the structure. The 
problem becomes thus substantially simpler and hence more economical. 

Since the motion the building starts from rest time that can 
selected origin 0), the initial conditions require that 


N 


Multiplying each Eqs. summing with respect and re- 


and 


Hence, whatever the values and the initial conditions for are 


Investigations California Special Publication No. 201, Dept. 
Commerce, Coast and Geodetic Survey, Washington, D.C. 


and 
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Substituting the value from Eq. 18b Eq. the velocity 
takes the form: 


Whenever and approximate solution the earthquake prob- 
lem obtained means the first modes, the ratios— 


rn (t) > [a; (20) 


quantitative idea the accuracy with which the initial conditions 
have been satisfied. (Although must vanish for physical reasons, 
sometimes convenient approximate means curves whose first 
derivative does not vanish the origin.) 


Whenever elementary analytic function, the solutions Eqs. 
satisfying the initial conditions Eqs. are given 


For example, 
Xosin wot 

gives 


and 


and the corresponding displacements become, Eqs. 


N 


and 
j#k 


Xo Ak [sin wo t + wo t cos wo t]; (wo = wr) . (24b) 


superposition, the displacements due half-sine pulse ground dis- 
placement— 
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—are given 


p(t) = 


The solutions Eqs. may used connection with each sine compon- 
ent the expansion into Fourier series and are practical import- 
ance connection with the graph maximum amplitude functions the 
frequency given Frank for each grade the modified Mercalli 
earthquake scale. 

When not elementary function, and particularly when 
known experimentally, the integration Eqs. conveniently performed 
numerical methods, shown Sect. analog machines. 


INFLUENCE INTERNAL DAMPING 


Although little known about the type internal damping developed 
actual structures, experimental evidence that damping steel occurs 
the coefficient viscous damping any given mode were proportional the 
frequency that mode. Since viscous damping easily handled mathemati- 
cally, this result assumed apply complete structure, and viscous term 
added Eqs. 16: 


The value the logarithmic decrement the jth mode— 


upper value (seldom attained) 0.9; and, hence, upper value for the 
coefficient viscous damping 


log, 0.9 


The solution Eqs. when sine function 


which 


Amplitudes and Periods Associated With Earthquake Preliminary 
Frank Neumann, Supplement, Quarterly Engineering Seismology {ST-75, October, November, 
December, 1950. 

and Cooper, Bell Telephone System Tech. Publications, Monograph 1561 (or Proceedings, Soc. for 
Experimental Stress Analysis, Vol. No. 1948.) 


0 
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a; wo (w?; = wo) 


Inasmuch the largest stresses due earthquakes occur, general, 
short time after the beginning the earthquake, and since the logarithmic 
decrement nearly the influence damping elastic earthquake stresses 
often negligible. Its essential feature damp out, rapidly, the higher 
modes vibration. 


was assumed the preceding sections that the motion the foundation 
was purely horizontal. When building rests elastic soil, instead, the 
earthquake displacements will produce, not only horizontal displacements 
the floors, but also rotation the foundation about axis through its centroid 
perpendicular the plane motion. (It assumed that relative horizontal 
displacement takes place between foundation and soil.) This additional 
motion usually called the the building. 

Calling the rotation the foundation radians, indicate (elastic 
stiffness coefficient for rocking motion) the reactive elastic moment the soil 
per unit depth the building when value was computed 
Biot’ terms the elastic constants E(Young’s modulus) and (Pois- 


son’s ratio) the soil: 


which the width the foundation. 

The modes the elastic structure rocking its foundation are conveni- 
ently determined means the concept systems” evolved 
Hans Bleich, “complementary system consider the 
building which the “shear springs” are frozen: System rigid beam 
elastically restrained against rotation system 
B,” consider the building which the “shear springs” are free act, but 
which rocking prevented. The modes system are, therefore, the 
modes, previously determined, the building rigid soil. 

The differential equation free vibration the rigid building (system 


which the frequency given 


and Experimental Methods Engineering Seismology,” Biot, Transactions, 
ASCE, Vol. 108, 1943, 365. 

Analysis Beam and Girder Floors,” Hans Bleich, Transactions, ASCE, Vol. 115, 

1950, p. 1023. 


and 


and 


EARTHQUAKE 


the moment inertia the structure about the axis rotation, indicat- 
ing the height the ith floor above the foundation (Fig. 2). The corresponding 
mode straight line 


N 


Center Line 


Axis Rotation 


Mr. has shown that the modes the rocking structure can 
represented linear combination and thus: 


The coefficients this linear combination are given 
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the value defined Eq. and letting 


N 


=1 


follows that 


and 


that must satisfy the condition 


the are normal modes. 
The frequencies the rocking structure are the roots the 


Once the values have been found from Eq. 43, the Eqs. and 
define and for each value and Eq. gives the new normal modes 
the rocking structure. (An analysis Eq. shows that lower than 
both and and that any other bound two successive 

The generalized co-ordinates the rocking structure satisfy the first 


second the equations 


and 


depending upon whether friction neglected taken into account; and the 
initial conditions are 


The absolute displacements the rocking structure are given, usual, 


now interesting notice that the stresses the rocking building are 
due exclusively the shear components the displacements z;. These are 
obtained subtracting from the total modes the rotational components 


| 
(45) 
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The shear displacements are then simply given 


When the number stories the building greater than the 
writing and the solution the determinantal equation for the frequencies 
(Eq. becomes lengthy task requiring particular care. simpler such 
cases obtain the frequencies and the modes the classical procedure 
successive approximations developed Vianello and Stodola, which 
makes use the fact that the modes are the deflections the elastic structure 
loaded its own inertia forces. 

determine the first mode the Vianello-Stodola 


(1) Assume deflection which possibly not too different from 
(2) Apply the structure the horizontal loads 


k 
(3) Determine the shear the kth floor: 


i=N 


(4) Determine the displacement the kth floor, relative the 1)th 
floor: (in which the spring constant under the 
floor) 

(5) Determine the new displacement the ith floor summing the 


k=1 


The displacement better approximation which can now 
used the initial deflection obtain nearer approximation and the 
process continued until the difference between two successive approximations 
becomes negligible. 

practice, since unknown, the structure loaded with loads 
and the deflections are times smaller than they should 
be. Hence the ratio any the corresponding gives 
approximate value 

The best possible value obtained taking the ratio the weighted 

N 


4 i=] 


t=! 


| | 
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The normalized mode given, terms the last approximation 
Eq. 

The second mode may obtained the same procedure, provided the 
initial approximation chosen orthogonal and the successive ap- 
proximations are made orthogonal easy see that, given 
any function the function— 


N 


used obtain functions orthogonal means the rth approxima- 
tion which, general, will not orthogonal 

The determination the third mode may similarly started means 
any function and the formula 


N N 


which makes orthogonal both and 

order able follow the convergence the method practical 
normalize the approximations each step. Since the results converge 
more rapidly when the initial approximations are near the true modes, pays 
obtain good initial guesses. These may obtained, for example, from the 


AND FREQUENCIES THE SHEAR BEAM 


Refer the building rigid soil vertical z-axis with the origin the 
foundation (Fig. 2), and let: 


the height the ith floor above the foundation; 
the height the ith story; 


the average mass per unit length (assumed constant) 
the average spring constant per unit length (assumed con- 
stant); 


the height the building; 


the mass the building; and 
the common value when does not vary from story story. 


Letting 
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take the form: 


and, dividing and taking the limit approaches zero, the partial 
differential equation so-called becomes 


this equation and (and hence their ratio k/m) are assumed 

independent (The solution the corresponding equation for linear 

variation k/m has been obtained the writer connection with further 

studies the dynamics shear When k/m constant, 
easy that the normal modes the shear beam are given 


and that the corresponding frequencies are given 


Merit ASCE, has proved, moreover, that the 
number floors larger than the first few frequencies, and, hence, the 
first few modes, the lumped building, for all practical purposes, are sufficiently 
well approximated those the uniform shear beam. 

Hence, when the ratios are approximately constant, one may 
obtain good approximations the lower frequencies and modes the building 
using for Eq. 55, its average value— 


—and for Eqs. and its average value 


obtaining 


Eighth International Cong. Theoretical and Applied Mechanics, Istanbul, Turkey, 1952. 
Pericds Uniform Cantilever Jacobsen, Transactions, ASCE, Vol. 104, 
402. 
433. 


and 
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Moreover, and are independently constant (uniform building 
constant story height), Eqs. and reduce to: 


shown the illustrative example Sect. these formulas give fairly good 


approximations the lower modes and frequencies building even 


not rigorously constant. 


THe NUMERICAL INTEGRATION THE EQUATIONS MOTION 


pointed out Sect. the equations motion (Eqs. 27) must 
integrated numerically whenever not simple analytical function. 
Numerous methods step-by-step integration are available for second-order 
equations, but few are simple and powerful those advanced 
and Fox and Goodwin,” and applicable Eqs. and 
27. 


integrate Eq. 


Eq. multiplied through the square the time spacing At, the 


4 
second derivative approximated the first two terms its 


central difference and both sides the equation are then operated 


2 


—Eq. reduces the simple recurrence equation 


= — Prart ¢ prt F,+ 


Fox 


12*Same New Methods for the Numerical Integration of (Sotoow Differential Equations,” 
and Goodwin, Proceedings, Cambridge Physical Soc., Vol. 45, Part 1949, 373. 
Difference Sheppard, Proceedings, London Mathematical Soc., Vol. 
31, 1899. 

actions, ASCE, Vol. 116, 1951, 590. 


and 
and 
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which the truncation error takes the form 


The truncation error the order and the convergence the process 
very rapid. 

convenient and accurate evaluate means the initial conditions ex- 
pressed terms averaged central differences. Neglecting differences 
order higher than the fourth, was done obtaining Eq. 63, Eqs. yield 


and 


and 


Solving Eqs. and Eq. 65a for and 


evaluated the required points, may obtained differencing the 
These values are introduced Eq. 63, each step, obtain improved valves 
seldom that this correction need used practical cases. 

Substituting for the derivatives the first term their central difference 
the following recurrence formula for the integration Eqs. 
similarly obtained: 


and 


EARTHQUAKE 


which 


and 


Eq. the coefficient viscous friction. The initial values and 
are both equal zero this case, and the truncation error proportional 
ILLUSTRATIVE EXAMPLE 

Modes and Frequencies Building Rigid Soil—The methods pre- 
sented the previous section will illustrated means the simple example 
five-story building, whose geometric and elastic characteristics are given 
Table 


i in. 107? in. ( 


The roots the corresponding determinantal equation, Eq. are given 
Table they were evaluated Newton’s method and synthetic division 
about one hour. Solving Eqs. successively with these values taking 
X;(5) the are and, Eq. the modes are found and 


entered Table 


Shear Beam 


2 
(1) (2) (3) (4) (6) (7) (8) (9) 
1.44 7.44 11.65 6.84 587 —587 12.17 
1.44 6.00 15.53 7.98 954 10.66 
1.44 4.56 15.53 9.12 1,101 587 12.18 
1.44 3.12 15.53 10.26 1,248 13.70 
1.68 1.68 15.53 11.40 1,395 20.72 
57.11 7.557 7.84 3.7 
424.89 20.61 23.53 14.1 
46.87 70.58 50.5 
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58, For SEVERAL VALUES 


Values of j equal to: 


—0.165 

0.041 
0.151 —0.025 —0.164 
0.144 0.16 0.101 


—0.166 0.149 
—0.107 
0.095 —0.116 
0.149 0.076 
0.102 0.130 


The uniformly distributed mass and spring constant the five-story build- 
ing are, Eqs. and 56, respectively, 0.0992 per in., and 
1.37 The frequencies the first five modes the correspond- 
ing shear beam, computed from Eq. 59, appear Table together with their 
percentage errors. They show that, although the ratios (Table 
are not constant, the first frequency fairly well approximated the funda- 
mental frequency the shear beam. The modes the shear beam, computed 
Eq. 58, appear for comparison Table 

Table which self-explanatory, shows the evaluation the first mode 
the Stodola-Vianello method, according Sect. starting from assumed 
deflection The second approximation, has errors not more 
than unit the third significant figure, when compared with the true mode 
Table 3(a). The corresponding frequency, evaluated Eq. 49, 
3,265; and 7.559—and has error less than 
0.03%. 

The evaluation the second mode the Stodola-Vianello method appears 
the initial approximation The first approximation then obtained 
the same operations used Table derive from The dis- 
example already orthogonal and hence (This indi- 
cates that, when the approximate deflection near the true mode, one need 
not orthogonalize the solution each step.) Repeating the process, 
found differ from the true Table 3(a) more than three units the 
third significant figure—that is, less than 3%. The corresponding value 


the frequency after the second approximation 


20.72—with error 0.53%; and, after the 


third approximation, 20.63, with error 0.05%. 


Floor 
1 2 3 4 5 
5 0.165 —0.093 0.020 165 0.165 
4 0.149 0.160 | —0.055 157 —O.151 
3 0.118 —0.120 0.116 135 0.114 
2 0.080 —0.039 | —0.168 101 —0.057 
1 0.040 0.127 0.139 057 —0.008 
| 
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TABLE THE FIRST 


0.070 0.400 10.26 0.0390 0.078 1.203 
0.039 0.439 11.40 0.0385 


TABLE THE SECOND 


(1) | (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) 
5 —10| —0.153 | —0.152| —1.767| —1.767 | —0.258| —0.338| —0.159| —0.159 | —1.852| —1.852]| —0.271 
4 0 0 0.001 0.018 | —1.749 | —0.219| —0.080|} —0.038| —0.038 | —0.582| —2.435| —0.305 
3 0 0 0.001 0.014) —1.734 | —0.190 0.139 0.066 0.066 1.019 | —1.416] —0.155 
2 10 0.153 0.154 2.386 0.652 0.064 0.330 0.155 0.155 2.409 0.993 0.097 
1 10 0.153 0.153) 2.381 3.033 0.266 0.266 0.125 0.125 1.944 2.937 0.258 ‘ 


The third mode was not computed since did not contribute substantially 
the solution, shown Sect. 9c. 

Modes and Frequencies Rocking soft elastic soil with 
10,000 per in., and Poisson’s ratio has elastic stiffness co- 
efficient for the rocking motion unit strip building wide, given 
32; thus: 1.145 The moment inertia the building, 
about axis through the centroid its foundation, Eq. 35, with and 
frequency the rigid building, Eq. 34, 7.043 49.60). The 
corresponding mode computed Eq. 36, appears Table together with 
the first two modes and the shear building rigid soil, taken from 
Table 3(a). 

The modes the rocking structure have been computed approximately 
means these three modes only. From Eq. and Table 0.998; and 
0.049; and with 57.11, 424.9 (from Table the frequency 


5 0.179 | 0.165 —0.166 0.173 0.077 0.162 0.129 
4 0.145 0.149 | —0.046 0.147 0.069 §£ 0.040 0.014 
3 0.110 0.118 0.095 0.114 0.055 » —0.099 —0.119 
0.075 0.080 0.149 0.077 0.037 
0.041 0.040 0.102 0.018 —0.102 —0.109 


Ki 
(2) (3) (4) (6) (7) (8) (9) 
11.65 0.105 0.105 6.84 0.0153 0.158 1.836 
15.53 0.124 0.229 7.98 0.0287 0.142 2.210 
15.53 6.5 0.101 0.330 9.12 0.0362 0.114 1.765 
a a = - = 4 
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vz: 
(16) 


(16) (17) (18) (19) (22) (23) 


equation, Eq. 43, becomes 


49.60 0.996 0.0024 


The first two roots this equation are 


0.998 


= 0.868 


0.049 


—and, Eq. 42: 49.60 57.11(0.868 424.9(— 0.003 26.57, 
from which 0.536; 0.465; and 0.002. Repeating the com- 
first two modes the rocking structure are given by: 


and 


(10) (12) (13) (14) (15) (1) 
0.290 1.836 0.268 2.764 7.640 89.0 0.164 
0.315 4.046 0.507 2.495 6.225 96.7 0.148 4 
0.201 5.811 0.637 1.988 3.952 61.4 0.118 3 
0.093 7.013 0.684 1.351 1.825 28.3 0.080 2 
0.023 7.611 0.668 0.668 0.446 6.9 0.040 1 
(14) (24) 
—0.377 1.655 5 
—0.106 0.175 4 
0.199 0.615 3 
0.354 1.950 2 
0.258 1.031 am 
and 
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The shear modes the rocking structure are evaluated setting 
and 

Motion Building Rigid Soil for Sinusoidal Foundation Displace- 
ment.—In order evaluate the influence the various modes the motion 
the building and gage the convergence the Noumerov-Fox step-by-step 
integration method, the building Table was subjected sinusoidal 
ground displacement unit amplitude and period 1.2 sec (the period 
important Fourier component actual earthquake record). With 
5.236, and 27.42, with and selected from Table and 
and from Table 3(a), Eq. gives the displacements: 


SINUSOIDAL FouNDATION DISPLACEMENT 


: sin wot sin wot 
(Mode (Modes (all sinwit | sin wet ri r2 rs 


(2) (3) (4) (6) (7) (8) (9) (10) 


2.204 2.173 2.178 0.123 —0.070| 0.000 
1.866 1.884 1.878 —0.070 0.062 —0.196 0.000 


Cols. and Table 7(a) contain, respectively, the coefficients 
the forced term sin when and modes are taken into account Eq. 
24. (These last coefficients have been obtained independently from rigorous 
solution the same problem the Laplace transform method.) These re- 
sults indicate that the present problem the contribution the higher modes 
the forced vibrations negligible and that permissible take into ac- 
count only two modes. Cols. and Table 7(a) contain the coefficients 
sin and sin the contributions the first two modes free vibration. 

was noticed Sect. that the initial conditions rest for the co-ordi- 
nates cannot rigorously satisfied unless all the modes are taken into ac- 
and shows that the rest conditions are very poorly satisfied the 
lower floors unless one takes into account all the modes. This expected, 
since the jump between and cannot well approximated 
Fourier series. The values which should identically zero, have been 
computed numerical check. noteworthy (and also lucky) that the 
solution this problem not too sensitive the initial conditions. 
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The Step-by-step Integration the Lagrangian displace- 
ments the structure due sinusoidal motion the ground were evaluated 
Sect. means the solution Eqs. given Eqs. 23. The values 
for appear Table for values between and 1.5 1.8 sec. 


TABLE INTEGRATION THE 
LAGRANGIAN FoR p(t) 


Approx- Approx- Approx- Approx- 
imate True imate True imate True imate 


(a) INTEGRATION 


1.335 


(b) INTEGRATION 


0.752 0.748 —1.118 
0.749 0.744 —1.214 
0.537 0.533 —1.192 
0.320 0.318 —1.059 —1.055 


0.054 0.054 —0.849 


—0.843 —0.839 0.076 
—0.823 —0.823 0.106 
—0.724 —0.726 0.162 


—0.767 —0.769 0.264 
—0.862 —0.862 0.427 
—0.991 —0.989 0.638 


—— 


The values these functions the same points have also been evaluated 
the Noumerov-Fox step-by-step procedure. For with 57.11 
0.312. With these constants and starting values, Eq. gives the 


values Table 8(a). comparison with the true values indicates satis- 
factory agreement although about one eighth the period 
Eq. 63, the results are These results agree satisfactorily 
with the true values after seventy-two steps, although the spacing the 
order the period 0.305 sec. 

The seventy-two steps the integration took approximately min., 


once the sine functions had been tabulated. The error was evaluated but 


True 
Pr pr Pr Pr Pr Pr Pr Pr 
(1) (2) (3) (1) (2) (3) (1) (2) (3) (1) (2) (3) 
0 0 0 5 —2.940 15 3.197 3.179 
1 0.032 0.048 6 —2.156 | 16 2.275 2.276 
2 0.312 0.335 7 —0.466 | 17 0.574 0.593 
3 0.884 0.901 8 1.480 18 —1.173 | —1.147 
4 1.507 1.508 9 2.885 | — _ _ 
0.004 0.006 19 55 0.870 0.866 
0.041 0.044 20 56 1.081 1.075 
0.135 0.138 21 57 1.236 1.229 
0.293 0.295 22 58 1.311 1.304 
0.505 0.505 23 59 1.296 1.291 
0.744 0.742 24 60 1.202 1.199 
0.974 0.970 25 61 1.056 1.056 
1.157 1.151 26 62 0.894 0.895 
1.265 1.258 27 63 0.750 0.753 
I 1.284 1.278 28 64 0.652 0.654 
1.221 1.217 29 65 0.610 0.610 
1.099 1.096 30 66 0.615 0.614 
0.951 0.950 31 67 0.645 0.641 
14 0.813 0.815 32 68 0.666 0.660 
15 0.716 0.718 33 69 0.643 0.637 
16 0.673 0.675 34 70 0.553 0.548 
17 0.680 0.680 35 71 0.387 0.383 
18 0.717 0.715 36 72 0.154 0.154 
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not used since did not influence the results within the accuracy the com- 
putations. 

Influence was stated Sect. that the influence 
damping the earthquake displacements consists, essentially, reducing the 
importance the higher modes. This revealed the case the structure 
under study, logarithmic decrement equal 0.9 assumed, leading 
value the damping coefficients (Eq. 29): 0.017 7.557 0.127; and 
0.017 20.61 0.350. Solving Eqs. with these values the 
damping coefficient f;, for sinusoidal ground motion unit amplitude and 
period 1.2 sec.—and using the modes and the undamped 
structure Table 3(a)—the floor displacements become 


sin wot 7.311 sin (wo 2°31’) 


comparison this equation with Eq. Sect. shows that the coeffi- 
cients the forced vibration are practically unchanged. The amplitude the 
higher modes, instead, rapidly damped out, and the frequencies are practically 
unchanged friction. 


10. 


The methods outlined this paper permit the solution dynamic problems 
shear buildings whatever the applied forces displacements may be. They 
must compared with other methods the basis efficiency and flexibility. 

direct measure their efficiency has been gained the writer solving 
the problem the structure without rocking three different ways: 


(1) direct integration the equations motion, using the Laplace 
transform; 

(2) numerical integration the same equations, using extension 
the Noumerov-Fox and 

(3) modes and numerical integration, using the methods this paper. 


The time involved the three solutions was roughly the ratio 3:3:1; but 
this ratio times does not indicate the main advantages the use modes— 
(a) the possibility simplifying the solution using limited number 
modes, and (b) the small amount additional work necessary analyze the 
structure under various types external loads displacements. 

The solution the rocking building direct numerical integration 
and requires min for each step, whereas the use comple- 
mentary systems permits the evaluation the rocking modes remarkably 
easy manner. The direct analytical integration the same problem, whenever 
possible, very cumbersome. 

The introduction internal damping does not present difficulties the 
present methods, but complicates the other two types solution considerably. 
seems possible state, therefore, that the methods presented herein may 


Comprehensive Numerical Method for the Analysis Resisting 


Whitney, Anderson, and Salvadori, Journal, Vol. 23, No. September, 1951. 
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more efficient and flexible than others used date. The availability such 
procedures should encourage the designer abandon the simplified and un- 
realistic methods analysis suggested codes (such the method) 
with obvious gains economy and safety. 

Finally, analogous methods are applicable tall buildings, which direct 
stresses are prevalent, and they may extended include the elasto-plastic 
behavior structures under dynamic loads. 
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